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Abstract 

Configurational states that are to be associated, according to Goldstein, with the basins in the potential energy landscape 
cannot be characterized by any particular basin identifier such as the basin minima, the lowest barrier, the most probable 
energy barrier, etc. since the basin free energy turns out to be independent of the energies of these identifiers. Thus, our 
analysis utilizes basin free energies to characterize configurational states. When the basin identifier energies are monotonic, we 
can express the equilibrium basin free energy as a function of an equilibrium basin identifier energy, as we explain, but it is not 
necessarily unique. 
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It is well known that most supercooled liquids (SCL) become viscous when their configurational entropy becomes 
negligible as they are cooled, provided the corresponding crystal is not allowed to nucleate. Our current understanding 
of glassy behavior is still far from complete, even after many decades of continuous investigation. In order to better 
understand the flow properties of viscous fluids, Goldstein proposed the potential energy landscape (PEL) picture 
using classical canonical ensemble 1, 2] to qualitatively discuss an interesting and sufficiently tractable scheme to 
study SCL and the glassy state by drawing attention to the potential energy basins that are expected to play a pivotal 
role in the thermodynamics of viscous fluids at low temperatures. The PE surface contains many local potential energy 
basins, each possessing a basin minimum (BM). The study paved the way to exploit PEL to study stationary SCL 
by the use of the partition function (PF) formalism. The stationary state is obtained under infinitely slow cooling of 
the disordered equilibrium liquid below the melting temperature Tm- The basins have other points of interest such 
as the lowest energy barrier (LB) or the most probable energy barriers, i.e. the barriers with the highest population 
(MPB) that most probably play an important role in crossover dynamics or diffusion, as the temperature is raised. 
We will collectively call these basin points of interest basin identifiers. Goldstein provided a qualitative description 
of the nature of the resulting landscape created by the distribution of basins and their minima, and drew attention 
to the idea of configurational states, which should be associated with individual basins and should change with the 
temperature T. Instead of being determined by the particular topology of the basin, the configurational state of 
the system associated with that particular basin, according to Goldstein, is uniquely specified only by the potential 
energy £ of the basin minimum. This will certainly make sense if the equilibrium value £ — £(T) of £ is in one-to-one 
correspondence with (i.e., is a monotonic function of) the temperature T. We wish to emphasize that throughout this 
work, we are only interested in non-negative temperatures. 

As we will see below, there is no thermodynamic requirement for £{T) to be monotonic. Even il£{T) is monotonic, 
is expressing thermodynamic quantities as a function of £{T) over a wide range of temperature an oversimplification, 
since it makes the shape and topology of the basins irrelevant? For example, if we consider a particular basin, then 
its free energy changes with T, but the basin minimum energy remains unaffected by T. Similarly, the free energy also 
changes as the system moves from one basin to another at the same temperature. The BM energy £ has no guarantee 
to remain the same. This change in £ is certainly not going to be reflected in £(T). Thus, we expect not only the 
equilibrium £(T) but all equilibrium quantities to have separate intrabasin and intcrbasin contributions. Indeed, this 
happens as we will see here [see (|31|) below]. We find that the equilibrium value of every thermodynamic quantity 
including £^ is a function not only of T, but also of the basin free energy. Only at low temperatures do we expect the 
latter dependendence on the free energy to become weak as assumed by Goldstein. In that case, the use of £ may 
lead to erroneous description if care is not exercised. What about other points like LB and MPB in the basins? Will 
they play a more important role at higher temperatures? Is there a thermodynamic basis for their importance? 

It is obvious that the shape and topology uniquely determine the free energy of the basins, which then uniquely 
determine the thermodynamics of the viscous fluid. Should we not describe configurational states by specifying the 
basin free energy? The latter is certainly monotonic in T. We note in this regard that, based on an analogy with spin 
glasses, the use of basin free energy in place of the BM energy, as originally advocated by Goldstein, has been proposed 



use of PEL has now become common in many disparate fields like glasses, proteins and clusters 9] , and has established 
itself as an important thermodynamic approach in theoretical physics. Thus, it is highly desirable to answer the above 
questions. It is this desire that has motivated this investigation. 

We briefly review the basics of the PF formulation in statistical mechanics jlO| as some of the steps in the evaluation 
of the PF are going to be crucial in our discussion here. The PF Z{T) is a sum of the Boltzmann weight exp(— /3i?„) 
for each configuration of energy En 



here n indexes the configurations of the system. However, as we can change the value of n associated with a configura- 
tion at will by relabeling it, the index n is not unique and, therefore, its possible values for equilibrium configurations 
are not very useful for thermodynamic description. Therefore, we need to relate n to a unique thermodynamic prop- 
erty of all configurations, and not just the equilibrium configurations. The property must be such that its equilibrium 
value is a monotonic function of T, so that it can be used in place of T if we so wish. It is not hard to identify such 
a property. Since each configuration has a unique energy, which is a thermodynamic quantity, we can use the energy 
En in place of n as the summation index. As the energy of interactions £"„ is coupled to the inverse temperature /3, 
see ([3): the equilibrium energy 
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is a monotonic function of T because of a non-negative heat capacity: E and T are in a one-to-one correspondence. 
Thus, the first step is to rewrite the PF in 

Z(T)^^t^(£;)e-'^^, (2) 

E 

where W{E) represents the number of configurations of energy E and defines the microcanonical entropy S{E) = In 
W{E). The sum in Z{T) at a fixed T is over all allowed energies E and, by construction, each summand is also 
a function of E. These energies and, therefore, S{E) do not depend on T. The second step in the evaluation 
is to recognize that the value of Z{T) for a macroscopic system, which is what we consider here, is determined 
by the dominant term in the sum. The dominant term corresponds to the equilibrium energy E at which the 
summand in ^ is maximum. Thus, Z{T) ^ W[E]e-l^'^. It is the presence of W{E) in Q that makes the summand 
strongly peaked at E. The determination of E for a macroscopic system is simplified by noting that E is almost 
a continuous variable for a macroscopic system. In terms of S{E), E is given by the location of the minimum of 
the free energy function F{T, E) — E ~ TS{E) at fixed T. In equilibrium, the canonical entropy S{T) = S{E) and 
free energy F{T) = F{T, E) = E{T) — TS{T) become functions only of T. The conditions for the minimum are 
[dF{T,E)/dE]T = 0, and [d'^F{T, E)/dE^]T > leading to 

[dSiE)/dE]^^^^^^ = (3, dS{T)ldT > 0, (3) 

which are always satisfied in equilibrium. 

At a given temperature T, only those configurations that have the energy E — E (or within a narrow width around 
it, depending on the heat capacity; we will neglect this width here) determine the thermodynamics through the 
entropy S{T). All energies other than E and, therefore, all configurations not included in W{E) are irrelevant at 
T. Thus, thermodynamics is highly selective. We will call configurations with E — E equilibrium configurations to 
distinguish them from all configurations. All thermodynamic quantities related to equilibrium configurations will be 
similarly called equilibrium quantities to distinguish them from quantities related to general configurations that may 
not be equilibrium ones. In the following, we will make this distinction carefully. 

As said above, the equilibrium values of E are strongly peaked about E; all these configurations have the same 
minimum free energy F{T). These configurations belong to different basins with different BM energies. Therefore, E 
at a given temperature must be independent of £. However, there is no guarantee that the the equilibrium values of 
£ are also strongly peaked about £. In other words, there is requirement that equilibrium configurations all belong to 
basins having their minima that are strongly peaked about £. It is certainly possible for equilibrium configurations 
to belong to different basins that have their minima at different energies, but still correspond to the same free energy 
F{T) because of the shape of the basins. In this case, it does not seem reasonable to classify equilibrium configurations 
only by the minima of the basins they belong to. We must also use the basin free energy in addition to the BM energies 
for this purpose as discussed above. However, it appears highly plausible that it can be done at low temperatures as 
originally proposed by Goldstein P, Q since the basin free energies are close to the BM energies in values. 

In his analysis, Goldstein has considered basin minima as the basin identifier, and has listed two conjectures that 
were common in the field 2] at the time: the basin PF zy,{T) is (i) independent of the basin's minimum energy £, 
and (ii) insensitive to the basins being explored. Utilizing these assumptions, Goldstein has expressed the PF as a 
product of the basin and BM PF's 

Z{T) = 2b(r)ZBM(r); (4) 

here Zh for a given basin is defined by considering shifted energies E — £ with respect to the minimum energy £ of 
that basin; see also Schulz . Goldstein has emphasized that basin anharmonicity or the curvature at its minimum 
|12|| may be very important. These are included in Zb, so that it is determined by the entire basin topology. According 
to Goldstein, all equilibrium basins have the same equilibrium basin free energy fh{T) = —Tin z^,. The BM-PF is 
defined as 

ZBuiT) = ^7VBM(f )e-'^^. (5) 

£ 

Here, iVBM(i^^) represents the number of basins whose BM are at energy £. The equilibrium BM energy £ = £{T) 
is the value of £ at which the summand in jSJl is maximum. The conditions for the maximum in terms of the BM 
entropy S'bmI^) = 1ii-/Vbm(i^^) are given by 

[dSBM{£)/d£]^^^ = f3, d£/dT > 0, (6) 
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which are the standard conditions of equihbrium; compare with (jSjl. It is clear that the BM description proposed by 
Goldstein ensures that £ is a monotonic increasing function of T. Since his approximation is expected to be good at 
low temperatures, we expect £ to monotonic increasing there. But we will see below that it need not be true at all 
temperatures, thereby limiting the usefulness of the BM-description to which we now turn. 

BM Picture. We take BM's as basin identifiers, but the discussion is easily extended to LB's or MPB's, or other 
basin identifiers. The PEL is a union of disjoint basins. A basin is indexed by j, and the lowest and highest basin 
energies are denoted by £j^ and £j, respectively, so that the basin does not exist outside the energy range AjE = 
{£j,£j). Since PEL does not depend on T, £j and £j do not depend on T. Let Wj{E) {E e ^jE) represent the 
number of distinct configurations of energy E in the j-th basin and introduce the entropy Sj (E) = In Wj (E) . It is 
obvious that Wj (E) is determined by the basin topology, which in turn is determined by the the interaction energies 
in the system. Let Wj{E) attain its maximum at Ej^ and introduce Aj^^E = {£j, Ej^). We now introduce the shifted 
PF 

z,iT)^ ^ W,{E)e-^(^~'^^ (7) 

EeAjE 

of the j'th basin. Let Ej{T), Sj{T) be the average basin energy and entropy, so that fj{T) = — Tlnzj(T) = Ej{T) — 
£j — TSj{T) represents the basin's shifted free energy. Of course, the conditions determining Ej(T) are 

[dS,iE)/dE]^^^^^^^ = /3, dE,iT)/dT > 0. (8a) 

Both conditions are always met. We wish to emphasize that Ej(T), Sj{T), and fj{T) do not represent equilibrium 
quantities yet; the latter are determined only after Z{T) is evaluated. (If each basin is treated as representing an 
independent system in a formal sense, then these quantities do represent equilibrium values for the particular basin.) 

We now group basins, indexed by j(A), into basin classes (BC) B\, indexed by A, so that all basins in a class have 
the same BM energy £ — £\. The basins in a class do not have to be close in the configuration space. The number of 
basins in B\ is Nbm{£\), and the corresponding BM entropy is Sbm{£\) = hi A'bm(£a)- The introduction oi Nbm{£x) 
and Sbm{£\) requires that they have a one-to-one relationship with £\. Let 

Za(T)= ^ zj(r), zx = Zx{T)INbm{Ex), (9) 

denote the shifted and the mean (per basin) shifted Ba-PF, respectively, so that 

Z{T) = e-^^>-Zx{T) = J2 e"'^^^+^™^^^^2A(T), (10a) 

A A 

£{T) ^ 5]£,e-^^^z,(T)/Z(r) ^ J] f,e-'^^^Z,(r)/Z(T). (10b) 
i A 

Here, £ = £{T) represents the equilibrium BM energy. The equilibrium free energy, entropy and energy are F{T) — 
-TlnZ(r), Sir) = -dF/dT and E{T) = F{T) + TS{T), respectively. 

At a given T, the probability Pr(r \B\) that the system will probe the BC Bx, i.e., any of the basins in B\ is 

Pr(T \Bx) = e-P^-Zx{T)/Z{T). (11) 

On the other hand, the probability that the system will explore a particular basin j among all the basins is 



Pr(rU)=e-'3^-z,(T)/Z(r), (12) 

whereas the conditional probability that the system explores the basin j, given that the system is confined to the BC 
Bx, is 

p,{T\\)=z,{T)lZx{T). (13) 

These probability distributions can be used directly to evaluate various entropies by using the Gibbs definition of the 
entropy as minus the mean logarithm of the probability distribution 

S = - <\nw>, (14) 
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where w is one of the probabihty distributions above. For example, the equilibrium entropy S{T) is obtained by 
considering the distribution in 112() : 

SiT) = -Y,P^iT\j)\nPTiT\j), (15) 

3 

The entropy S{T) due to the number of different equilibrium BC's, and the entropy S\(T) due to basins within a BC 
are similarly given by 

S{T) = - ^ Pr(T \Bx)\nPr(T\Bx), (16) 

A 

Sx{T)^-J2p,iT\X)lnp,{T\X). (17) 

J6J(A) 

It is easy to see that the temperature coefhcient d£/dT is given by 

T'^dE/dT = A£j{T)AEj{T) Pr(T \j) , (18) 
i=i 

where A£j{T) = £j—£, AEj{T) = Ej{T)—E{T), and A''b is the number of basins. The derivative is a cross-correlation 
between two fluctuations, A£j(T) and AEj{T,V). Since cross-correlations do not usually have a fixed sign, there is 
no theoretical reason for £ to be a monotonic increasing function of T. This does not mean that the temperature 
coefficient cannot be positive for many physical systems in a certain temperature range. The same argument also 
works for other equilibrium identifier energies. 

The sum over A in (|10allUb|l is over different basin classes, and the summand is also uniquely determined by each 
class. However, the index A is not a unique labeling. Therefore, we need to associate with each class a unique 
thermodynamic class property. We also need to make sure that its equilibrium value is monotonic in T. Thus, our 
task is to find a basin quantity that is not only a unique property of each class, but also determines the summand 
in HlUallUb() . One of the simplest choice is to take the basin minima £ as the basin quantity for the simple reason 
that each BC has a uniquely BM energy £\, and it also appears in the summand in HlUallOb|) provided we take z\ to 
depend explicitly on £\ in addition to T : = zx{£\,T). This contradicts the assumption by Goldstein 0. However, 
this is a common choice; see jl^ . If the interest is to study high temperatures near the dynamic crossover, one might 
take the LB-energy £j'^ of the j-th basin as the quantity of interest. In this case, one would group basins not by the 
their BM energies, by their LB energies, so that all basins in the BC A have their LB at £x^. Similarly, one can take 
the MPB-energy or some other special energy of the j-th basin as the quantity of interest. The analysis in all 

these cases would be carried out in the same manner as the analysis we carry out here for the BM picture, and would 
not be pursued further except for a few remarks later. 

Common Assumption. The assumption z\ = Z\{£\,T) as a function of two different variables T, and £\ makes 
the summand an explicit function of £\. Note that £\ in the summand in HlUallUb() is a temperature-independent 
quantity. Let us follow the consequence of this assumption. Let the mean basin free energy resulting from the mean 
basin PF z\ for the basin class B\, and the mean basin entropy be f{£\,T) = —T\nz\{£\,T), and S{£\,T) = 
-[df{£x,T)/dT]g^, respectively. The mean basin free energy f{£x,T) = ExiT) - £x - T S{£x,T) is obviously 
obtained by minimizing the free energy function fx{£x, E,T) = E — £x — T Sx{£x, E) with respect to E at constant 
T and £x- Here, Sx{£x,E) = Sx{E) — S^y\_{£x), and exp[S'A(i?)] represents the number of configurations of energy 
E that belong to Bx- The condition for the minimum is \dSx{£\-,E)ldE\^^ = P a.t E = Ex{T), which determines 
fi£x,T) = Ex{T)-£x-T S{£x,T)- here, S{£x,T) = Sx[£x, ExiT)]. 

The sum over the index A is now replaced by the BM variable £ — £x. Because of the assumed iS-dependence, 
the general summand in IjlUal lUb|) becomes an explicit function of £, and we can minimize the corresponding free 
energy function Fb{£, T) = £ + f{£, T) — TSbm{£) with respect to £ at fixed T to determine Z{T) for a macroscopic 
system. The minimum of Fb{£,T) is located by the condition [dFB{£,T)/d£]T — 0. This condition is satisfied at the 
equilibrium BM-energy £= £{T) = £{T), see (|10b|) . It is also equivalently given by the solution of 

[dSBM{£)/d£],^j = /3{1 + [idfi£, T)/d£)T],^j}, (19) 

and determines the equilibrium free energy F{T) = Fb{£,T) and the BM-entropy Sbm(T) = Sbm{£)- The solution 
of (|19|l yields £= £{T) as a function of a single variable T. The equilibrium mean basin free energy and entropy are 
fh{T) — f{£,T) and Sh{T) = S{£,T), respectively. It is easy to see that the form of the equilibrium free energy 
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Fb{T) = f{£,T) + £{T) — TSbm{£) is the same as the free energy obtained by Goldstein in except that the 
equations determining the equihbrium BM-energy are different; compare (jHl and 119(1 . The two conditions become 
identical if / is taken to be independent of £, as was assumed by Goldstein. 

Let us assume that £ = £{T) is a monotonic function of T, so that we can invert it to express T = T{£), which 
can be used to express all T-dependent quantities like fh{T), F{T), Sbm(T), etc. as functions of the equilibrium 
BM-energy £. We can also express f{£\,T) as a function of £\ and £. All this means is that T and £ are equivalent. 
The equilibrium free energy fh{£) = fh\T{£)\ is a single- variable function of £ and says nothing uniquely about the 
f -dependence in the two- variable function f{£,T). This is most clearly seen in the Goldstein's scenario. While fh{T) 
in his analysis is independent of £, we can still express fh{T) as fh{£)- To conclude that /b(f ) imphes that fh{T) is 
a two variable function having an additional iS-dependence would be a wrong conclusion. A similar inversion can be 
carried out in terms of the average LB-energy £j^^ (T) so that all quantities that depend on T can be expressed as a 
function of f (T) , which suggests that there is nothing unique in the choice of BM energies as the thermodynamic 
quantity to replace the sum over A. 

In the above analysis, we have only considered one {[dFB{£,T)/d£]T = 0) of the two conditions needed for 
the free energy minimization. The other condition is that the curvature of the free energy at the minimum be 
positive:[d^ Fb{£ ,T) / d£^]T > 0. This issue has been investigated elsewhere and will not be pursued here. 

Current Analysis. We now proceed to prove our central result that £\ (or any other basin identifiers) cannot be 
the choice of the thermodynamic quantity to replace the sum over A for two different reasons that will be discussed 
below. We observe that W{E) is a sum over various basins containing the energy E: W{E) = J2j-E£A e ^j(^)- We 
also observe that even for basins in a single BC B^, Zj{T) need not be equal. 

1. Evaluating Zj{T). We first prove one of the assumptions by Goldstein |3| that Zj cannot depend explicitly 
on the BM energy £j, see though it most certainly depends on the shape of the basin, i.e. on j. For example, the 
basin curvature and not £j determines the vibrational frequencies (including anharmonicity) which, in turn, determine 
the free energy fj{T) = — Tlnzj(T) in the harmonic approximation. The latter is measured with respect to £j, so it 
must be independent of £j . 

To be convinced that this is most certainly correct, consider two basins j —1 and 2 that are identical in all respects 
(shape and topology), but have different BM energies, which we assume satisfy £i < £2- Let us introduce SjE = E — £j. 
Let Wj{6jE), Sj{SjE) — lnWj{SjE) denote the number of , and entropy due to configurations of energy at a height 
SjE above the basin minimum £j in the j-th basin. It is obvious that 

W;{S,E) = W,{E). (20) 

Now, for diE = 52E = 5E, it is clear that S[{SE) — S2{SE) for the two identical basins, from which it immediately 
follows that dS[{SE)/dSE = dS2{SE)/dSE = (3. The free energies of the two basins are given in terms of the same 
average energy height 5E{T) = E.j{T) — £j, where Ej{T) is the average basin energy in the j-th basin, and the same 
entropy S[[SE{T)] = 5^[(5£;(T)] : 

MT) = SE{T) - TSr[5E{T)] = E,{T) - £, - TS,{T), 

and are equal, as expected, after we use (|20|) . Here, Sj{T) is the canonical basin entropy in the j-th basin. We can 
also check the veracity of the free energy equality by directly evaluating the shifted PF in jT)) . We will not do this as 
it is trivial. The extension to many topologically identical basins is trivial. Thus, we conclude that the free energies 
fj (T) do not depend on the BM energies £j . 

Indeed, a very general but direct proof can be given that the basin PF Zj{T) cannot depend on its BM energy £j. 
Thus, the corresponding basin entropy Sj{T) is also independent of £j. For this, let us shift all energies E ^ E = 
E — C hy some constant C in 0. The number Wj{E) of states, all having the same energy E, remains unchanged 
under the shift by C: Wj{E) Wj{E ) = Wj{E); compare with (^01) . Thus, Zj transforms under the shift as 

Zj{T) Y.E' W]{E)e-f^^^' see ©, and remains unchanged for any arbitrary C, including C = £j,£f, or £f^^ 
corresponding to the j-th basin. 

This central result requires adopting a different approach than the above for the PF-evaluation, which we present 
below. 

For a macroscopic system, Zj is determined by the maximum summand in |7J corresponding io E ^ Ej = Ej (T) G 
AjE. In terms of the free energy function (pj{E,T) = E — TSj{E), the energy Ej is determined by {d(pj /dE)T ~ 0. 
As usual, this condition yields 

{dS,iE)/dE)^^^ = /?, dEj{T)/dT >0,E-E A,„,i?; (21) 

the restriction on the allowed values of Ej ensures that we are only considering non-negative f3. Over this temperature 
range, Sj{T) = Sj{Ej) and Ej{T) are monotonic increasing. On the other hand, the shifted free energy of the basin 
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is /j(r) = Ej{T) - £j - TSj{T) is monotonic decreasing since Sj{T) = -[dfj{T)/dT]e^ > 0. Since there is no 
fj-dependence in fj{T), as we have argued above, there is no need to show that the derivative dfj{T)/dT is at fixed 

The energy landscape is topologicaUy very complex, with various basins very different from each other, even if they 
have their minima at the same energy. Thus, Ej{T) in different basins at the same temperature T will be usually 
different not only from each other, but also from E{T) that appears in (j^J. In this sense, all basins are independent 
at this step of the analysis as they all do not represent equilibrium configurations. By the last step of the analysis, 
only equilibrium basins will survive for which Ej{T) — E(T). 

2. Evaluating Zx{T). As shown above, Zj{T) is not a function of £\. We now need to replace the summation 
over j in @ by a suitable thermodynamic quantity. There are various choices for the sum in like the basin free 
energy fj{T), the basin energy Ej(T), or the basin entropy Sj{T), all equally good. All these quantities are monotonic 
in T. We will choose fj(T) in the following to be the thermodynamic quantity to replace the sum over the index j. 
As the basins in do not all have the same free energy, we classify each basin in (jSJ according to its free energy 
/ at a given temperature T. Let N\{f,T) denote the number of basins of free energy / at a given T in B\, and let 
Sx{f,T) = In Nx{f,T) be the corresponding entropy. For a macroscopic system at a given fixed T, Z\ is dominated 
by the basins in B\ for which F\(f,T) = f{T) — TS\{f,T) is minimum as a function of / at fixed T, the conditions 
for which are {dF\/df)T = 0, and {d'^F\/df'^)T > 0. The resulting entropy and free energy at the minimum (/ = f^^) 
are denoted by S\{T) — Sx{fx,T), and F\{T) = f\~ TSx, respectively. One of the conditions for the minimmn is 
given by 

[{dSx/df)T ]f^j^ = /3. (22) 

From ((T^ . we find the conditional probability Pj{T |A) = exp[— /3(/j — ^^a)]- Using this in ((TTIl . we can evaluate the 
entropy due to the number of basins in B\. It is easy to see that this entropy is precisely S\{T), as expected. 

Note that £\ is independent of T. Using once more the shift argument E E = E — C, and the invariance of Zj 
so that the sum in ^ contains the same number of basins N\{f,T), we similarly conclude that Z\{T) also remains 
unchanged for any arbitrary C, including C = E\. Thus, the shifted BC free energy Fx{T) is independent of £x, 
although it most certainly depend on the BC B\. In other words, it cannot be written as a two-variable function 
F\{£x,T). The corresponding entropy S\{T) is also independent of £\. 

Even for the singular case, when Zj{T) for each basin in the BC Bx is the same, all basins [their number being 
Nbm{£\)] contribute. In this case, Nx{f,T) — Nbm{Sx) and the above formulation continues to work. 

3. Evaluating Z{T). Using the evaluated Zx in Z, we find 

Z = 5]e-'^[^-+^^(^)l. (23) 

A 

Since we have already seen that Fx{T) is independent of f a, the summation over A cannot be replaced by a summation 
over £x. We need to look for a BC quantity that can be uniquely related to the index A and which also uniquely 
determines the summand. For this we proceed in a standard manner as follows. Let Af{J-, T) denote the number of 
BC's at a given T with the same T = £x + Fx{T), where T is the unshifted BC free energy now with respect to 
the zero of E. Since F for different basin classes are now measured from the same common point = 0, it is clear 
that if we consider all those BC's that have the same free energy F at some T, then this F no longer depends on 
the individual BM energy £x. The PF Z is dominated by the BC's for which F — TS is minimum over F; here 
S{F^ T) = \iiJ\f{F^ T). One of the conditions for this minimum at F is 

[{dS/dF)^]^^y = P, (24) 

as expected. Hence, we finally conclude that the free energy of the system is given by 

F{T) = ~T\nZ = F -TS; (25) 

here, S{T) = S{F, T). It should be obvious at this point that the dominant contribution in (|23|l will usually mix BC's 
with different £x- 

The equilibrium free energy is the sum of three terms: F{T) — fhiT) ~ TSbc(T) — TS{T), where Sbc denotes the 
equilibrium basin class entropy and ^hiT) the equilibrium basin free energy. Expressing (phiT) — E\^{T) — TSh{T), 
we immediately see that 

S{T) = S^{T) + Sbc{T)+S{T). (26) 
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Introducing T\= £\ + F\(T), we can express the probability Pr(T \B\) , see that the system will explore B\ 
in terms of J-x, 

PiiT\Bx)^e-^^^''^l (27) 

It is easy to check that using the above probability in (|16() . we obtain the above entropy S{T). In terms of the above 
probability in H27|l . the equilibrium BM energy, see (jlObp . is expressed as follows: 

£iT) = ^£xPr{T\Bx). (28) 

A 

Thus, unless S{T) = 0, so that there is only one BC, the equilibrium BM energy £{T) is not given by the BM energy 
of a single BC. In other word, the system will usually explore many basin classes of different BM energies, unless 
S{T) = 0. it is expected that S{T) vanishes at low temperatures. Thus, the use of BM's makes perfect sense at low 
temperatures as Goldstein had argued At higher temperatures, many BC's with different £\ are involved in 

(El- 

The last prediction is different from the BM approach in which the equilibrium basins all have the same average 
BM energy £{T). Let us try to understand this point carefully. Let us group the BC's according to their free energy 
JF. The number of BC's in a group G{J', T) corresponding to a given T is Af{J-, T). Let 

£{T,T)^ J2 £x/^i:F,T) (29) 
Aee(JP,T) 

denote the mean BM group energy in the group Q{J-,T). Then, 

f(T)= J2 ^^(^,T)e^(^'^)-'3[^-^l. (30) 

For a macroscopic system, the dominant term in the above sum H30|) corresponds to J-^ —F^ so that 

£{T)^£(T,T). (31) 

Unless N{J-^T) = 1, the sum in H29() contains many BC's, each with different £\. Since the summand in (|29|l does 
not have any coeeficient similar to W{E) in (jSJ, there is no reason for £'^s to peak around the mean £{T,T). Thus, 
the equilibrium BM energy £{T) represents a mean group BM energy over many different BM energies £\. 

From (|31|l . we conclude that £{T) should be properly thought of as the equilibrium value of a two- variable function 
£{T,T). In this form, the change in £{T) has two different contributions. One contribution is due to the intrabasin 
temperature change dT at constant J^, and the other contribution comes from basin changes caused by free energy 
change in dJ-' at constant T. There is no reason for either contribution to vanish in general. 

It is also easy to evaluate equilibrium values £^{T), or £^^^(T) using our formalism. Using the probability Pr(r |j) , 
see l)12|l . which is given by 

Pr(T|j) = e-W.(r)+^^-^m]^ (32) 
we can calculate the above equilibrium energies as follows: 

j j 

Following a similar analysis that wc have carried out here using BM as the point of interest, we can derive a similar 
relation like H31|l for the above two equilibrium energies. We will not do this here. By comparing the equilibrium 
energy E{T) with the equilibrium values in (|33|1 . we can draw conclusions about the importance of various basin 
identifiers at a given temperature. 

The landscape analysis developed here does not involve the equilibrium BM energy £, because £j are not the 
natural energies that appear in the PF. Despite this, it is possible to express thermodynamic quantities in terms of £ 
under certain conditions that we elucidate below. The equilibrium basins at a given temperature form a subset of all 
basins and determine the equilibrium thermodynamics. If £(T) is monotonic in T, we can invert the relation £(T) as 
before: T = T{£). Thus, we can express any function of T as a function of £. For example, the equilibrium basin free 
energy fh{T) can be expressed as a function ol £ = £{T). This is one way to express the equilibrium free energy as a 
function of £{T). There are other ways to do this, as we will discuss below. 
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We can use E{T) at a given T to select only the equilibrium basins whose average energy is E{T). This is precisely 
what is done when equilibrium configurations are generated in simulations carried out at a fixed T without any 
consideration of the basins From these equilibrium configurations, equilibrium basins can be identified. These 
basins also appear in (|28|l and determine £. But simulations can also provide information about whether there are 
many different BC's that contribute to £. (Unfortunately, to the best of our knowledge, simulations carried on 
finite-size systems have not been used to unequivocably answer whether the spread in the allowed values of £ are 
thermodynamically significant are not.) The equilibrium free energy f{T) of these basins will depend on the average 
topology of the basins. For example, in the harmonic approximation, /(T) will be determined by the curvature (given 
by the set i' of vibrational frequencies) at the BM, and should be expressed as a function J{i',T). These equilibrium 
frequencies will change with T. Being a function of T, they can also be expressed as a function of £ as discussed above. 
Thus, the equilibrium basin free energy can be expressed as a function /(5, T). This provides us with another way to 
express /(T) as a function of £. Such a dependence on £ of the equilibrium basin free energy should not be confused 
with any explicit fj-dependence of Zj{T) in |(7|), or Z\{T) in (|10a|l . since such a dependence has been shown not to 
exist. The free energy in H19(l is not this equilibrium free energy in its new form. Moreover, simulations usually do 
not generate basins other than the equilibrium basins. Thus, they provide no information about the free energy / 
that appears in (fT^ . 

It should be clear that depending on the complexity of the basin potential, there are many ways to express f{T) as 
a function of £. For example, the anharmonicity, to be represented in short by a set of parameters a, of the equilibrium 
basin will also depend on T. Now, we may express both u and a as a function oi £, or decide only to express v as a, 
function of f , and leave a as a function of T. This will give us two different forms of the equilibrium basin free energy, 
with both carrying the same information as a function of T, but not as a function of £. Thus, such representations 
are obviously not unique. This is because of the independence of the basin free energy and £. On the other hand, if 
a quantity has an explicit f-dependence, we can immediately obtain its equilibrium value at a given T by replacing 
£ by £. For example, S^m{£) can be expressed in a unique way as a function of T by identifying Sbm{T) = S-bm{£)- 

The above argument is easily extended to expressing thermodynamic quantities in terms of £^^ {T) , £^^^ [T] , or 
any other particular average energy related to the basins. There is nothing unique about the BM energy in the 
landscape, except possibly at low temperatures where the equilibrium free energy is close to £(T). 

In summary, we have made a distinction between basins and equilibrium basins and quantities related to them, which 
is usually not done. With this distinction in place, we have shown that the shifted basin free energy is independent of 
particular basin energies like that of BM, LB, MPB, etc. It is a function only of T. Thus, the equilibrium basin free 
energy is also a function of T, which is consistent with the first assumption of Goldstein. It is possible to express the 
equilibrium basin energy as a function of £ (T) ,£^^ (T) ,£^^^ (T) etc. in presumably many ways, as discussed above, 
but only if £ {T) ,£^^ (T) ,£^^^ (T) etc. are monotonic functions of T. The unshiftcd cquilibr ium BC free energy T 
that contributes to H25I) is also independent of the BM energy £ by construction. This is consistent with the second 
assumption of Goldstein, albeit in a modified form. Thus, the configurational states that are to be associated with the 
entire topology of the basins cannot be characterized simply by their BM energies £, LB energies £^^, MPB energy 
£;MPB gQ]-|2g other particular basin energy ; however, the basin free energy can be used to specify the configurational 
states. The current analysis deals only with the basin free energies. We have found that there is not one average BC 
of a given average BM energy £{T), but many different BC's are explored unless S{T) — 0. 

It is our pleasure to thank Andrea Corsi for his comments on the work. 



[1] M. Goldstein, J. Chem. Phys. 51, 3728 (1969). 

[2] M. Goldstein, J. Chem. Phys. 64, 4767 (1976). 

[3] M. Mezard and G. Parisi, Phys. Rev. Lett. 82, 747 (1999). 

[4] B. Coluzzi, G. Parisi, and P. Verrocchio, Phys. Rev. Lett. 84, 306 (2000). 

[5] B. Coluzzi, A. Crisanti, E. Marinari, F. Ritort, and A. Rocco, Eur. Phys. J. B 32, 495 (2003). 
[6] T. Aste, and A. Coniglio, J. Phys.: Condens. Matter 15, S803 (2003). 

[7] P.W. Anderson in "III Condensed Matter", Lecture Notes of the Les Houches Summer School, 1978, ed. by R. Balian, R. 

Maynard, and G. Toulouse (North-Holland, Amsterdam, 1979). 
[8] C. Dasgupta and O.T. Vails, Phys. Rev. E 59, 3123 (1999). 

[9] D. Wales, "Energy Landscapes", Cambridge, Cambridge University Press (2003). 
[10] L.D. Landau, E.M. Lifshitz, "Statistical Physics", Part I, 3rd Edition, Pergamon Press, Oxford (1986). 
[11] M. Schuiz, Phys. Rev. B 57, 11319 (1998). 
[12] C.A. Angell and K.J. Rao, J. Chem. Phys. 57, 470 (1972). 
[13] F.H. Stillinger and T.A. W eber, P hys. Rev. A25, 978 (1982). 
[14] P.D. Gujrati, |cond-mat704l2548l . 



9 



